Spatially Modulated Interaction Induced Bound States and Scattering Resonances 
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We study the two-body problem with a spatially modulated interaction potential using a two- 
channel model, in which the inter-channel coupling is provided by an optical standing wave and its 
strength modulates periodically in space. As the modulation amplitudes increases, there will appear 
a sequence of bound states. Part of them will cause divergence of the effective scattering length, 
defined through the phase shift in the asymptotic behavior of scattering states. We also discuss how 
the local scattering length, defined through short-range behavior of scattering states, modulates 
spatially in different regimes. These results provide a theoretical guideline for new control technique 
in cold atom toolbox, in particular, for alkali-earth- (like) atoms where the inelastic loss is small. 



Feshbach resonances (FR) and optical lattices (OL) 
are two major techniques in cold atom toolbox. FR can 
be used to control the interactions by tuning a bound 
state in the so-called "closed channel" to the scatter- 
ing threshold via magnetic field, laser field or external 
confinement At resonance, the s-wave scattering 

length diverges and the system becomes a strongly inter- 
acting one. OL can strongly modify the single particle 
spectrum of atoms, which suppress the kinetic energy so 
that the interaction effects are enhanced. With these two 
methods, many interesting many-body physics, such as 
BEC-BCS crossover, superfluid to Mott insulator transi- 
tion and strongly correlated quantum fluids in low dimen- 
sions, have been studied extensively in cold atom systems 
during the last decade Q. 

In this letter we theoretically study a new control tool 
for cold atom system. It is analogous to OL because it 
also makes use of two counter propagating laser fields 
that lead to a periodic modulation of laser intensity in 
space; however, its main effect is not acted on single par- 
ticle, but manifested on the interaction term when two 
particles collide with each other. It generates a spatial 
modulation of two-body interaction, i.e. the two-body 
interaction potential not only depends on the relative 
coordinate of two particles under collision, but also de- 
pends on their center-of-mass coordinate. As far as we 
know, this is a situation not encountered in interacting 
systems studied before, ranging from high-energy and nu- 
clear physics to condensed matter systems. One can ex- 
pect a spatial modulated interaction potential will result 
in many fascinating phenomena. 

The explicit model under consideration is schemati- 
cally shown in Fig. [TJ The open and closed channels 
are different orbital states and are coupled by laser field. 
Such an optical FR has been studied before for uniform 
laser intensity ^ , and has been observed experimentally 
for both alkaline and alkaline-earth- like Yb atoms ^] . In 
this letter we shall consider the situation that the laser 
field is a standing wave whose intensity, and therefore 
the coupling strength between open and closed channel, 
is modulated periodically in space. Such a setup allows 
one to control spatial modulation of intcr-atomic inter- 
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FIG. 1: A schematic of the model we studied, 
detail description. 



See text for 



action on the scale of sub-micron. Recently, this has 
been realized in ^^^Yb condensation 0, although the op- 
tical standing wave is a pulsed one. Alkiline-earth(-like) 
atom like Yb is particular suitable for such an experiment 
because the narrow ^Sq-^Pi inter-combination transition 
line can avoid large inelastic scattering loss. In this ex- 
periment, a spatially modulated mean- field energy has 
been observed from diffraction pattern in a time-of-flight 
imagine Q. However, the theoretical study of this sys- 
tem is still very limited, and even the two-body problem 
has not been studied. In this work we show that surprises 
indeed arise even in the two-body problem of this model. 

Coupled Two-channel Model. We consider a two-body 
Hamiltonian for a FR in which the open and closed chan- 
nels are modeled by two square well potentials Q 
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where R = (ri + r2)/2 and r = ri — r2. For r < tq. 
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In this model Vq is given by the background scattering 
length flbg as tan(fcof'o)/(fco?'o) = 1 — o-hg/fo where fco = 



2 
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FIG. 2: First four bound states as a function of the amplitude 
of coupling SI. The parameters for this plot are q = 0, Sc = 
0.05£r, JsTabg = -0.01 and Kro = 10"^ En = h^K^/m is 
taken as energy unit. 




FIG. 3: rw{x,r) (where w{x,r) is the "wannier wave func- 
tion") for the first four bound states, a — 2tt/K is the "lat- 
tice spacing". The parameters for this plot are ec = 0.05-Er, 
Kabg = -O.Ol, Kro = 10"^ and Kr = 0.1. 



\/mVo/h, and Vc is determined by the binding energy of 
closed channel molecule ec through Vc — H'^tt'^/ {mrQ) — ec- 
The size of inter-atomic potential ro is much smaller than 
all the other length scales. Conventionally, the inter- 
channel coupling f2(R) is a constant independent of R. 
Such a model captures all key features of a FR P, 0] ■ A 
bound state appears at threshold and causes scattering 
resonance at {Iq = i/ec/l/?], and the scattering length is 
given by 
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where /3 — 32roabg/(97r^). 

Now consider the situation V, depends on R. Solving 
the Schrodinger equation follows two steps: (i) in the 
regime r < tq, because the —fP'S/^/m term commutes 
with the Hamiltonian, we consider the wave function of 
following form 



sin kr , sin fcr , _ , 

i^o = a R ; Vc = KR) 

r r 

where a(R) and &(R) satisfy a coupled equation 
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a(R) + 0(R)6(R) = ea(R) 
6(R) + r2(R)a(R) = eb(R) 



(5) 

(6) 
(7) 



where e ^ E — h'^k'^ /m. There will be a set of eigen- 
function a;(R), 6;(R) and fc; that give rise to the same 
energy E. The eigen wave function in the regime r < tq 
should be assumed as 



, X ■^-^ , sin kir 

R,r =y"A — - 
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(ii) The superposition coefficient Ai, the binding energy 
E for bound states, as well as the phase shift d{E) for 
scattering states, are determined by matching the wave 
function in the regime of r > tq at r = tq for any R. 



Hereafter we will consider an explicit situation where 
ri(R) = il.cos{Kx) {x denotes the x-component of R). 
Note that there is still a discrete translation symmetry 
X — >■ x+2tt/ K, we can introduce a good quantum number 
"crystal momentum" q. In the regime r > rg, ipc — 0, 
and for the bound states whose energy E < fi^ / (Am) , 
'0^(x,r) can always be expanded as 



„-r,^{q+nK)^~4mE/h'^ 



Eq. © can be viewed as the Bloch wave function for 
molecules. And for the low energy scattering state whose 
energy is greater than but close to ft^5^/(4m), we have 



sin(A:r — S) 
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where k = ^ mE/h? — 9^/4, and (5 is a function of k. 

Results 1- Bound States: In contrast to the uniform 
case where there is only one bound state when f2 > f2o, 
in this case we find a sequence of bound states as f2 in- 
creases, as shown in Fig. [5J This is because the periodic 
structure of coupling r2(R) leads to a " band structure " 
for the molecules, and as the coupling strength increases, 
the molecules with zero crystal momentum but in differ- 
ent bands touch the scattering threshold one after the 
other. We can introduce the "wannier" wave function as 



w{x — xq, r) = 



K/2 



K/2 



(11) 



As shown in Fig. [31 the "wannier" function for the bound 
states that appear at larger fl have more oscillation, 
which means that they come from higher bands. This 
can also be illustrated from the symmetry of Un in the 
Bloch function of Eq. as summarized in the Table U 
for the first four bound states. The first two bound state 
has even parity while the other two have odd parity. 
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FIG. 4: The effective scattering length defined as Eq. ([T2|) 
tSoff/|t[bg| as a function of n/fio- (b) and (c) are enlarged plot 
around Q/ilo = 2.64 (b), 9.20 (c). The arrows indicate the 
positions at which we plot the local scattering length aioc in 
Fig. EJa-d). 



FIG. 5: The local scattering length aioc as a function of posi- 
tion x/a for r^/fio = 0.71, 2.55, 2.64 and 2.7 (a-d). The solid 
blue line is calculated results, the black dashed line is the fit- 
ting formula Eq (|17P or p8[) . and the green dash-dotted line 
in (a) is from simple replacement formula Eq. p5[) . 
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TABLE I: Symmetry of Bloch wave function for the first four 
bound states 



Results 2 - Effective Scattering Length: For the scat- 
tering state wave function, at large r only the first term 
in Eq. (|10l) will not exponentially decay, and the asymp- 
totic behavior of the scattering wave function is still the 
same as that in the uniform case. Hence we can introduce 
an effective scattering length as 



lim 

fc->0 



tan(5(fc) 
k ■ 



(12) 



Note that though the interaction is spatially dependent, 
the effective scattering length defined as Eq. is a 

spatial independent one. Among the first four bound 
states, floff only diverges when the second bound state 
appears at threshold, as one can see by comparing Fig. 
m^a) with Fig. This is because the divergence of 
implies the first term in Eq. ^TU\\ goes like l/r, which 
should be smoothly connected to a zero-energy bound 
state with non-zero Uq. Therefore, for the other three 
bound states whose Uq — 0, their coupling to the low- 
energy scattering states vanish and will not cause diver- 
gency of aeff. In Fig. UJ^c) we show that UoB diverges 
when the sixth bound state (whose Uq ^ 0) appears at 
scattering threshold, but the width of resonance becomes 
narrower compared to Fig. Hl^b) because this bound state 
comes from higher band and its coupling to low-energy 
scattering state ( i.e. the absolute value of Uq) is smaller. 
Results 3 - Local Scattering Length: At short distance 



the wave function Eq. (fTU)) satisfies the Bethe-Peierls 
contact condition and display l/r — l/aioc(2^) behavior, 
hence we can introduce a local scattering length as 



aioc (2^) = - lim 



ril)o{x,r) 



r^ro dr{r^-!o{x, r)) 



(13) 



Unlike in the uniform case, acff and aioc are different. 
Similar situation has also been encountered for scatter- 
ing in confined geometry [3], lattices ^ and mixed di- 
mension [9|. What is unique here is that oioc is spatially 
dependent. Naively, one may think that aioc(a;) can be 
obtained by replacing Q, in Eq. (j4]) by local ri(x), i.e. 



aioc (a;) = ahg 

~ abg [1 



PVL^ cos^{Kx) 
' ec + PVL^ cos^iKx) 
/3n^ cos'^{Kx)/ec] 



(14) 
(15) 



where the second line is valid for small fl. This formula 
in fact corresponds to an oversimplified approximation 
in our model that the kinetic energy term of the center- 
of-mass motion (— ?i^V^/(4to)) is completely ignored in 
Eq. (HI). In fact, what we really obtained from the wave 
function Eq. pil| is 
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1 - 2U2Cos{2Kx)/Uo 



- 2U2Kcos{2Kx)/Uo 
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The second line is also valid when ft is not too large, so 
the coefficient Um>2 is small enough that can be ignored. 

Away from a resonance, Kucs <^ 1, Eq. (jl6p can be 
well approximated as 



aioc (a;) = acff 
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In fact, we show in Fig. [Ha), (b) and (d) that the for- 
mula Eq. pT)) (dashed black line) is a very good approx- 
imation to the actual results (solid blue line). In Fig. 
[5]Ja) we show the simple replacement formula Eq. 
already significantly deviates from the actual results in 
weak coupling regime. From Fig. ^h) and (d) one can 
also see that the mean value of aioc(a;) changes sign as 
Ceff changes sign. At resonance, a^^ — > 0, Eq. (fT6|) can 
be approximated as 



aioc(a;) = — 
K 



1 - 



Un 



2U2 cos{2Kx) 



(18) 



We show in Fig. [SJc) that Eq. (IT5|) is also a very good ap- 
proximation to actual aioc at resonance. Hence, we show 
that oioc behaves very differently in the regime nearby or 
away from a scattering resonance. 

Implications to Many-body Physics: In summary, we 
have revealed a number of novel features in the two-body 
problem with a spatially modulated interaction potential, 
which have strong implications for many-body physics 
and provide new insights for developing new tools for 
quantum control in cold atom systems. 

First, when flcff diverges, the system enters a strongly 
interacting regime and is expected to exhibit univer- 
sal behavior, which can even be manifested in the high 
temperature regime [loj . For a two-component Fermi 
gas, it provides a new route toward BEC-BCS crossover 
physics, and "high-temperature" superfluid may exist in 
this regime. The periodic structure will add new ingre- 
dient to the crossover physics. 

Secondly, for the low-energy states whose energy ^ 
E^, the energy dependence of scattering length can be 
ignored and the many-body system can be effectively de- 
scribed by a pseudo-potential model: 
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2m ^ 
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where R^j — (r^ -I- rj)/2 and r^. 



(19) 

■j. It is very 

important that aioc(R) in the pseudo-potential of Eq. 
(fT9|) is given by Eq. ([T6|) from the two-body calculation, 
so that a two-body problem of the Hamiltonian Eq. (fT9|) 
can produce correct low-energy eigen-wave function and 
the effective scattering length as from model potential. 

For bosons, with a mean-field approximation, Eq. (jl9p 
implies that the interaction energy should take the form 
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aioc{x)n^ {x)dx 



(20) 



which leads to a modulation of condensate density n{x) 
and self-trapping nearby the minimum of aioc(a^)- It is 
very likely a strong enough modulation of condensate 
density will eventually result in the loss of superfluid- 
ity and the system enters an insulating phase. If so, it 
provides a completely different mechanism for superfluid 



to insulator transition where the transition is not driven 
by suppression of kinetic energy as in conventional OL. 

Final Comments: In this work we choose a coupled two 
square-well model whose advantage is that the physics 
can be demonstrated in a simple and transparent way. 
However, some more sophisticated effects in real system, 
such as the inelastic loss, are ignored. We have also im- 
plemented more systematic scattering theory which in- 
cludes these effects and found that the physics discussed 
here will remain qualitatively unchanged. These results 
will be published elsewhere 11|. 

Moreover, the formalism used in this work can be eas- 
ily generalized to other realizations of spatial modula- 
tion of interactions. For instance, in a magnetic FR, one 
can consider the presence of a magnetic field gradient so 
that the closed channel molecular energy varies spatially. 
This effect is particularly important for a narrow reso- 
nance. One can also optically couple the closed channel 
molecule to another molecular state via a bound-bound 
transition, which leads to a periodic variation of molecule 
energy jl^]. Similar effects as discussed in Results 1-3 



also present in these cases |11| . 
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Appendix: In this appendix, we present some details 
of solving the two-body Schrodinger equation. Using the 
discrete translation symmetry, we expand 



a„ and &„ satisfy coupled matrix equation 



4m 



^n+l) ~ 
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This matrix has a set of eigen-values and their eigen- 
vectors fe/'„}- Hence there are a set of wave func- 
tions sharing the same energy E, which ar 



_ ^iqx Sill(fc,^r) ^ iriKxa 
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where — yJm[E — ef) /h is a function of E. In general, 
the eigen-states take the form 



c,l 
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\ 'fin / 



where 



<„(r)=^A^ 



sin(fc;'r 



For bound states, to match the boundary condition 
with Eq. ^ at tq in both open and closed channels, we 
obtain a matrix equation M^i{E)Al = where 



ML,i = (fc? cos(fc/ro) + \j{q + nK) 



AmE 

2- -^sin(/c«ro))a^„ 



Therefore for a given q the eigen-energy E is determined 
by Det(Af) = and 



/ h^(q + nK) 



^A,aisin(fcVo) (28) 



For the scattering states, ; and M|„ ; (n ^ 0) are 

the same as bound state, while 



Mq ; = sin(fc; ro) cos(A;ro — 6)k — cos(fc; ro) sin(fcro — 5)aQ 



where k — ^JraEjW- — cp- j^. In this case Det(Af) = 
gives rise to the relation between phase shift b and energy 
E. Uj[ (n 7^ 0) is also the same as Eq. (l28l) . while for 
n = 0, 
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